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Abstract: We address the question of defining the second quantised monopole creation 
operator in the 3+1 dimensional Georgi-Glashow model, and calculating its expectation 
value in the confining phase. Our calculation is performed directly in the continuum 
theory within the framework of perturbation theory. We find that, although it is possible 
to define the “coherent state” operator M(x) that creates the Coulomb magnetic field, the 
dependence of this operator on the Dirac string does not disappear even in the nonabelian 
theory. This is due to the presence of the charged fields (W^). We also set up the 
calculation of the expectation value of this operator in the confining phase and show that it 
is not singular along the Dirac string. We find that in the leading order of the perturbation 
theory the VEV vanishes as a power of the volume of the system. This is in accordance with 
our naive expectation. We expect that nonperturbative effects will introduce an effective 
infrared cutoff on the calculation making the VEV finite. 
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1. Introduction 

Common lore has it that the condensation of magnetic monopoles is responsible for confine¬ 
ment in nonabelian gauge theories [1]. Numerous lattice investigations of the monopoles 
in pure Yang-Mills theories are available in the literature [2], Nevertheless, there are many 
open questions regarding the status of magnetic monopoles in confining theories. The 
main problem with the monopole condensation scenario is that no gauge invariant defi¬ 
nition of the monopole operator can be given in most theories of interest. Moreover, no 
gauge invariant observable corresponding to magnetic charge exists at least in the SU(2) 
pure Yang-Mills theory [3]. Thus serious doubts remain regarding the significance of the 
monopoles for confinement in pure gluodynamics. 

On the other hand there undoubtedly exist theories where magnetic monopoles can be 
given a proper gauge invariant meaning. A prime example of such a theory is the Georgi- 
Glashow model in 3+1 dimensions. The model comprises of the SU(2) gauge field coupled 
to the adjoint Higgs 

L = -\f 1 2 u + (D ix H) 2 -V(H 2 ). (1.1) 

The theory is believed to have two phases as a function of the Higgs potential V(H 2 ). In 
the Higgs phase, where the VEV of the Higgs field is nonvanishing, the magnetic monopoles 
exist in the spectrum as heavy particles. Closer to the boundary of the confining phase 
they become light, and it is natural to expect that they condense in the confining phase of 
the theory. The magnetic charge (which, as opposed to the pure Yang Mills theory, does 
exist here as a physical observable) is then expected to be spontaneously broken in the 
confining phase. 

This behaviour has been all but established in the supersymmetric cousin of the Georgi- 
Glashow model [4], It is however worth noting, that the confinement in the deformed N=2 
supersymmetric theory is essentially abelian, and thus differs from the confinement in QCD 
in many important aspects [5]. Confinement in the Georgi-Glashow model on the other 
hand is expected to be fully nonabelian, and thus bear close similarity to QCD. 
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Since the magnetic charge certainly exists as a physical observable in the Hilbert space 
of the Georgi-Glashow model, one is tempted to try to use this model as crutches on the way 
to the pure Yang-Mills theory. More specifically, one would like to start with the Georgi- 
Glashow model and trace the fate of the monopoles when the mass of the Higgs field in the 
confining phase is taken to be very large. In this limit the theory becomes the pure SU(2) 
Yang-Mills theory. Of course, the perturbation theory, which is valid in the Higgs phase, is 
not valid deep in the confining regime anymore. Nevertheless, one can hope that the heavy 
Higgs field can be integrated out perturbatively so that one can at least understand what 
becomes of the monopole creation operator in the pure Yang-Mills limit. Since we know that 
the magnetic charge disappears in this limit, two options seem to be open for the monopole 
creation operator. The first one is that it simply reduces to some run-of-the-mill operator 
with the vacuum quantum numbers. This operator could conceivably be a convenient 
order parameter for confinement, even though it does not probe symmetry breaking any 
more. The other possibility is that the condensate of the Georgi-Glashow monopoles deep 
in the confining phase is so stiff, that it costs the energy of the order of the Higgs mass to 
excite it. In this case any excitation coupled to the condensate will be very heavy, and in 
the limit the monopole creation operator would simply reduce to a unit operator. In this 
case, it would become useless for discussions of the properties of pure Yang-Mills theory. 
In particular it could not be an order parameter for the deconfining phase transition at 
finite temperature. The former option seems more likely, since generically one expects light 
scalar excitations (glueballs) to couple to the “radial” part of the magnetic monopole field, 
and thus to excite the condensate. Nevertheless the second option is a logical possibility 
worth exploring/ruling out. 

To study this question one has to start from the beginning, and the beginning appears 
to be the construction of the monopole creation operator in the continuum Georgi-Glashow 
model. Surprisingly, although there exist a large number of lattice simulations in the pure 
Yang-Mills theory, there have been to our knowledge no attempts to explore this aspect of 
the Georgi-Glashow model. The continuum literature is also sorely missing on this point. 
There have been attempts to construct the monopole operator in the Georgi-Glashow 
model [6], but they appear to be incomplete. The aim of the present paper is to construct 
a gauge invariant monopole creation operator in the Georgi-Glashow model, and to set up 
the perturbative calculation of its expectation value in the confining phase. We construct a 
continuum gauge invariant operator which creates the Coulombic magnetic field (without 
the Dirac string) of the monopole and has the correct commutation relation with the 
magnetic charge density: 

[M(x),B k (y)\ = - j- — M(x ), 

g (,x - y) 6 

A t yr 

[M(x),j^{y)\ = — 5 3 (x - y) M(x). (1.2) 

This operator is constructed by analogy with lower dimensional theories as a “coherent 
state” creation operator, that is an exponential of an operator linear in electric fields. An 
important property of the monopole operator is that it is not rotationally invariant even 
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though the magnetic field created by it is spherically symmetric. This is directly related 
to the fact that the Dirac string is present in the definition of the monopole operator. The 
direction of the string has to be specified and it affects the commutators of M with local 
gauge invariant quantities, like the spacial components of electric current. We stress that 
these commutators are not hopelessly divergent at the position of the string, but merely 
depend on its orientation. 

We are also able to set up the perturbative calculation of the expectation value of M 
in the confining phase. We find that, as in the calculation of similar averages in lower 
dimensions [7], the relevant path integral is dominated by a classical configuration. This 
configuration interpolates in time between two states with magnetic charge one half (in 
units of the magnetic charge of the ’tHooft-Polyakov monopole). At time plus infinity the 
magnetic charge density is spread out homogeneously over space, while at minus infinity 
the state contains a point-like ’tHooft-Polyakov monopole and a balancing charge minus 
one half spread out homogeneously in space. Viewed from a four dimensional perspective, 
the configuration looks like radial magnetic current emanating from the endpoint of a stem 
(world line of the ’tHooft-Polyakov monopole). We therefore call this configuration “a 
dandelion”. The physical amplitude described by this configuration corresponds to the off 
diagonal matrix element between two components of the vacuum wave function (Ti|M|To), 
where |'f'i) is a state with magnetic charge one, and |To) a state with magnetic charge zero. 
This is a typical situation in theories with a spontaneous broken symmetry (in the present 
case the symmetry generator being the magnetic charge). A surprising property of the 
leading dandelion configuration is that it is not spherically symmetric as one might naively 
expect. This is the direct consequence of the lack of rotational invariance in the definition of 
the monopole operator. The problem of string dependence is not new and is well known in 
the framework of the lattice gauge theory [8]. A possible solution to it has been suggested 
by Frohlich and Marchetti [9]. Although the prescription of [9] directly applies to the 
calculation of the correlation functions, and not to the operator itself, roughly speaking 
it corresponds to the definition of M as a weighted average over operators which create 
two Z 2 strings rather than a Coulomb magnetic field. This prescription goes outside the 
framework of the coherent state operator, which we have adopted here and is much more 
complicated. We therefore do not consider it in this paper, and we believe that it will not 
alter qualitatively our results. 

Although the structure of the leading configuration suggests spontaneous breaking of 
magnetic charge, we still find (M) = 0 even in the confining phase. The reason for this is 
that the classical action of the dandelion is infrared divergent. In particular, as a function of 
the infrared cutoff (the size of the system V ) we find (M) oc exp{—^ In V}, where c is a pure 
number of order unity. This is not unexpected for the simple reason that in the perturbation 
theory the gluons are massless. In the presence of massless particles the theory has no 
scale, and thus we expect a logarithmic divergence of the classical action. The situation 
in the perturbation theory is thus similar to the Kosterlitz-Thouless type phases in 1+1 
dimensional theories [10], where continuous symmetries are “almost” spontaneously broken 
and correlators of the order parameter decay not exponentially, but rather as a power of the 
distance. Beyond perturbation theory in the Georgi-Glashow model we expect the situation 
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to change. Nonperturbative effects will introduce an infrared scale, the confinement radius, 
and we expect this scale to cut off the infrared divergence of the classical action and to 
lead to a finite VEV of M. 

This paper is organised as follows. In Section 2 we discuss the construction of the 
monopole operator. In Section 3 we set up the calculation of its expectation value and 
discuss the properties of the dandelion configuration. Section 4 contains some brief com¬ 
ments. 

2. The monopole creation operator 

The classical held configuration of the ’tHooft-Polyakov monopole is 

r° f (A r b 

H a = h{r )— , H“ = ^e aib -, (2.1) 

r g r A 

where h( 0) = /(0) = 0, lim,.-^ h(r) = v, lim^oo /(r) = 1 and v is the vacuum expectation 
value of the Higgs held. The distance scale on which the helds approach their asymptotic 
values is given by their masses. Since we are only interested in the infrared properties of 
the theory, we will look for an operator which creates a point-like monopole and thus from 
this point on hx h(r) = v, f(r) = 1. 

This configuration carries one unit of magnetic charge. The gauge invariant magnetic 
charge in the Georgi-Glashow model can be dehned in the following way [1]. First one 
dehnes a gauge invariant held strength as 

VM = - i Hj3 £»■” H‘ (D ll H) b (D„H) C , (2.2) 

in terms of the nonabelian held strength tensor 

F^ = d li Ai-d v A^ + ge abc A b li Al, 

and covariant derivatives of the Higgs held H a 

(D^H) a = d^H a +ge abc A b fl H c . 

The ’tHooft tensor can also be written as 

= d(H a A a u ^-d v (H a T“) - i e abc H a d^H b d u H c , (2.3) 

with the unit vector held H a = H a /\H\. The magnetic current is then dehned as 

kf?(x) = ±e^ Xa d" F^(x). (2.4) 

This dehnition is convenient since the magnetic charge density so dehned for all smooth 
configurations is equal to the topological charge density of the Higgs held. The total 
magnetic charge is then equal to the topological charge of H. 

Qm = ^ J d 3 x C (x) = J d 2 Si e ijk e abc H a d j H b d k H c , (2.5) 
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where the last integral is over the surface surrounding the zeros of the Higgs field. It 


should be noted that the definition eq. (2.4) is not unique. One can add to this current 
any other conserved current with the same quantum numbers whose charge vanishes on 
the ’tHooft-Polyakov monopole configuration. In particular it is useful to consider 


jjf = e^\ad v 


H a (x)F Xaa (x) 


( 2 . 6 ) 


On the ’tHooft-Polyakov monopole configuration the charge corresponding to this current 
is the same as Qm since the covariant derivative of the Higgs held vanishes outside the 


monopole core. On the other hand the density of this charge, unlike that of eq. (2.5), is not 
constrained to vanish outside the points where H = 0, and thus behaves like the majority 
of other local operators. For this reason we prefer to use this alternative definition of the 
magnetic current in this paper. 

Our aim is to construct a quantum operator that creates a monopole in the fully second 
quantised theory. This operator has to satisfy the following commutation relations 

47T 


[M(x), jo (y)] = — 5' 6 (x - y) M{x ) 


(2.7) 


It is natural to constrain the operator M further by requiring a somewhat more restrictive 
commutator 


M(x)B k (y)M (x)t = B k (y) + - 


i {x - y)k 


( 2 . 8 ) 


9 (x- y ) 2 ’ 

where B k = \e kimH a Fj 1 m . This latter commutator ensures that the operator M creates a 
Coulomb like magnetic field. 

These commutation relations are reminiscent of the ones satisfied by soliton creation 
operators in 2 and 3 dimensions [7]. There is however one crucial difference, and that is 
that the monopole is a nonlocal object, and thus has a nonlocal commutation relation with 
local gauge invariant fields, e.g. H,(x). 

The most natural way of going about the construction of M is to try a coherent state 
like operator, that is an exponential of the linear functional of electric fields. We thus start 
with the following ansatz 

M(x) = D(x)M A (x ), (2.9) 

where 

M a {x) = exp (i [ d 3 y\i(x - y)H a {y)E^(y)\ , (2.10) 


( 2 . 11 ) 


with A i the classical vector potential of a point-like Dirac monopole: 

1 4 

A i(x) = - eij -4- (cos 9 - 1) . 

9 n 

Here we have chosen the Dirac string to run in the direction of the negative X 3 axis, and 
have defined = ( xi,X 2 )■ The operator M A creates a vector potential of the Dirac 
monopole which in the isospace is oriented in the direction of the Higgs field. In particular 

(x - y)i 


M A {x)Bi(y)M ] A (x) = Bi(y ) + 


(x - y ) 2 


- 47 TZi 6 2 (x ± - y±)0(-(x 3 - y 3)) 


( 2 . 12 ) 


- 5 - 








Note that the operator Ma is fully gauge invariant. It does however create the magnetic 
field of the Dirac monopole including the contribution of the Dirac string. As a result it 
does not create a magnetic charge, since the magnetic flux due to the Coulomb magnetic 
field is balanced by the flux of the Dirac string. The role of the operator D in eq. (^1]) is to 
create the magnetic held of the negative Dirac string alone, and thus cancel the unwanted 
singular magnetic held in eq. (2.12). 

To construct the string operator D we turn for inspiration to the classical monopole 
solution eq. (2T). This solution is written in the gauge where the vector potential is regular 
everywhere except at the location of the monopole. One can however transform this into 
unitary gauge, where the Higgs held is constant everywhere in space. The transformation 
that achieves this is given by 


U(x) = exp{—i- a ■ </>} , 


(2.13) 


where 9 and (j) are azimuthal and polar angles respectively, cr are Pauli matrices and the 
unit vector 


b l = |x_l|<9 i4> = -a 


The formal gauge transformation with this gauge matrix gives 

W(x)Ai(x)U(x) + — U\x)dj,U (x) = \i(x)a 3 , 
9 

W(x)H(x)U(x) = (73, 


(2.14) 


(2.15) 


where we have defined Ai = A^a a and H = H a cr a . The transformation eq. ( |2.15 ) is not a 
proper gauge transformation. The rotation matrix U is singular along the negative X3 axis, 
where the direction of rotation <f> is undefined. As a result the transformed configuration is 


not equivalent to eq. (2.1), as the magnetic held corresponding to eq. (2.15) has an extra 


Dirac string, and the magnetic charge of eq. (2T) is cancelled by the string contribution. 


The second quantised version of the transformation eq. ( 2.13 ) (or rather its inverse) is 
precisely what we need to dehne to be able to get rid of the Dirac string contribution in 
eq. (|2.9|). It is however somewhat cumbersome to work directly with the transformation 


eq. ( 2.13 ) as it changes the gauge variant variables everywhere in space. It is more con¬ 
venient to redefine it by a proper gauge transformation which achieves the same result as 
eq. ( |2.13 ) everywhere except infinitesimally close to the Dirac string. 

Dehne 


Q 

U{x) = exp{—i- a ■ $} , 


(2.16) 


where the function 9 is equal to the azimuthal angle 9 everywhere except inside an infinitely 
narrow tube around the negative X3 axis, e.g. 


9 = 9 


! - - sign(x 3 ))/(x_i_) 


/(0) = 1, /(x > e) = 0. 


(2.17) 
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The action of this transformation on the configuration eq. (2.1) gives 


where 


H = HU = 03 cos(# — 9) + a ■ r± sin(# — 9 ), 

Ai = U\x)AAx)U(x) + — UUx)diU(x) 

9 

= — - ^5j(/>[cos(0 — 9) — cos 9]H — di{9 — 9)P\ — di(f)sm{9 — 0)P2 ) , 


(2.18) 


P 1 = a ■ cj), P r 2 = cos (9 — 9) a ■ r± — sin(0 — 9) <73 . (2-19) 

In the formal limit 6 = 9 the configuration eq. ( |2.18 ) becomes equivalent to eq. ( 2. 15| ) , 


however at any nonzero value of the regulator e in eq. ( 2.17 ) it is fully gauge equivalent to 
eq. ( 0 ). 

The transformation 
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V(x) = U\x)U{x) = exp{— i —-— a ■ 0} 


( 2 . 20 ) 


is gauge equivalent to U(x), but acts even on gauge variant fields only in the infinitesimal 
neighbourhood of the negative X 3 axis. To construct the operator that creates the Dirac 
string we must understand how to write down the second quantised version of the transfor¬ 
mation eq. ( 2.20| ). The difficulty here is in the fact that in writing down the matrix V we 
have assumed explicitly that the field H on which it is acting points in the third direction 
in the isospace. For any other background, V should be additionally rotated by a regular, 
field dependent transformation. It is however not very transparent how to deal with a 
transformation whose parameters themselves depend on quantum fields. The alternative is 
to take the limit e —> 0, in which case the Higgs field does not change at all under the action 
of the monopole operator. The vector potential however is affected by the transformation 
V even in this limit. Clearly, if H does not change, the vector potential classically must 
acquire a non singlevalued piece if the resulting magnetic field is to have a nonvanishing 
divergence. It is not difficult to write down a multivalued vector potential which gives a 
magnetic field of a Dirac string. 


i-A = t4 = 0 , 


= - 4>5{r)9(-X3) H a , 
9 


( 2 . 21 ) 


where the radial delta function is one dimensional rather than two dimensional to preserve 
the correct dimensionality of the vector potential. A shift of A by /x is affected by the 
following second quantised operator 


D(x) = exp 


d 3 y (f>{x-y) 5(|x ± -y±|)q-^- 0(-(x 3 -y 3 )) H a (y)E°-(y) ) . (2.22) 


|X_L - y_L 


Even though this expression looks suspect, it is in fact a well defined operator in a non- 
abelian theory. On the microscopic level (ultraviolet cutoff scale) the electric field operator 
in a nonabelian theory has a discrete spectrum which is quantised in units of g. Thus the 
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apparent ambiguity in the definition of cj) is not felt by any physical observable since it only 
leads to ambiguity of an integer multiple of 2tt in the phase of eq. ( 2.22| ). The operator 
D(x ) is gauge invariant and is similar in form to M A (x). In the following therefore we 
will not write it down explicitly. We will continue denoting the vector potential of the 


monopole by A i but will understand it as the sum of Aj of eq. ( 2.11 ) and yi of eq. ( 2.21 ). 
Hence, with this definition the monopole creation operator satisfies the correct commu¬ 


tation relation with the abelian magnetic field, eq. (2.8). This commutator is independent 
of the direction of the Dirac string, even though the definition of M uses explicitly the 
Dirac string in a particular direction. However the dependence on the string does not dis¬ 
appear from all commutators of M with gauge invariant operators. In particular consider 
the following commutator 

M\{x) [(£? {y)f - B?(y)\ M A (x) = [(B?(y)) 2 - B 2 {y)\ + 2e ijk Bf (y)Xj(y - x)(D k H(y)) a 


+ 


~ x)(D k H(y))° 


(2.23) 


This explicitly depends on the vector potential A*. 

We note that the commutator eq. ( |2.23| ) has a direct analog in the abelian theory. 
Consider an abelian gauge theory with a charged scalar field 0. The action of an abelian 
monopole creation operator on the spatial components of electric current cjfDi<f> is 

M\x) [pDifty) - Di</>*(i>(y )] M(x) = [<j>* D^y) - D^* <j>{y)\ + i\(x - yW{y)^{y) . 

(2.24) 

Thus the action of the monopole creation operator produces a long range electric current. 
Similarly the gauge coupled kinetic term of the scalar field transforms as 

M\x)Di4>*Di<j){y)M{x) = D^*D^y) - i\(x - y) ((/>* Di<p(y) - D^* </>(y)) + A 1$*c\>. 

(2.25) 

It is this property of the magnetic monopole operator in theories with charged matter fields 
that leads to dependence of the operator on the choice of the string, as has been discussed 
extensively in [8,9]. 

The commutator in eq. ( 2.23| ) is very similar in physical meaning. The operator DiH 
is naturally identified with the charged vector gauge boson fields W± (the two components 
of the vector field perpendicular to the Higgs). The square of the nonabelian magnetic 
field contains the minimal coupling of the charged vector field to the abelian gauge 
potential. Thus eq. (2.23) is simply a transcription of eq. (2.25) in terms of spin one 


charged matter. The potentially disturbing property of the commutator eq. ( 2.23 ) is that 
it is much more nonlocal than the commutator with the abelian magnetic field. Whereas 
the abelian magnetic field created by the monopole operator decreases as l/(x — y) 2 , the 
nonabelian magnetic field decreases only as l/|x — y\. As we shall see in the next section 
this nonlocality has an important effect on the calculation of the expectation value of the 
monopole creation operator in the confining phase. 


3. Perturbative calculation of (M) 

In this section we set up the calculation of the expectation value of the monopole ere- 


















ation operator in the confining phase. Although we do not expect the perturbative result 
to be reliable, we do expect to see substantial differences between the behaviour of the 
expectation value in the confining and Higgs phases. 

In the Higgs phase, since the monopoles are massive particles, the correlation function 
of the monopole creation operator at large distance should behave as 


(M*{x)M{y)) oc exp ( - y,\x - y \), 


(3.1) 


where y is the monopole mass. Thus the expectation value in the finite but large volume 


is 


(M) oc exp ( — fiL). (3-2) 

In the confining phase we expect the operator M to have a finite expectation value. 
However, since confinement is intrinsically nonperturbative, we still expect to find a van¬ 
ishing expectation value in the perturbative setup. However we do expect the volume 
dependence at finite volume to be much milder. 

Consider the path integral calculation of (M): 


(M) = / DHDA exp ( - S[A\) M[A\ 


DA ex p 


- \ / d 4 x(F; u -f; u f 


(3.3) 


with 


foi = \i(x)H a 5(x 4 ), fij = 0. (3.4) 

We are interested in calculating the expectation value deep in the confining phase, where 
the Higgs field is very heavy. In the first approximation we therefore neglect the Higgs 
contribution to the action. The c-number field eq. ( |3.4|) is the shifted Dirac vector potential, 
which multiplies the operator of the electric field in the exponent eq. (2.1C). The extra 
c-number contribution to the action f 2 arises when writing the expectation value in the 
functional integral form as discussed in detail for example in [11]. 

Path integrals of this type arise frequently in calculating expectation value of opera¬ 
tors which create topological excitations [7]. At weak coupling the path integral can be 
calculated in the steepest descent approximation. This becomes obvious once one realises 
that a simple rescaling of fields A —* jjA, leads to the appearance of the factor 1/g 2 in 
front of the action in eq. (|3.3| ) . 

Our aim is therefore to find a classical configuration which minimises the action in 
eq. ([37^). First we cast the path integral expression into a more intuitively appealing form 
by changing variables 

A ^( x ) -» A l( x ) - X ti (x)9(-x 4 )H a , (3.5) 

with Ao = 0 and the spatial components of A* given by the sum of eq. ( PHD and eq. (|T2l| ). 

Under this transformation 


flD 


- d lvK]H a + A [u {D^H) a 9(-x 4 ) - A[„ S /l]4 H a S(x 4 ) 


The path integral for the calculation of VEV of the monopole operator becomes 
(M) = [ DA exp j -- fc - - (a u D^H - X^hY 0(-® 4 


(3.6) 


(3.7) 
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(3.8) 


where f,, v is the magnetic field of a point-like ’tHooft-Polyakov monopole at X 4 < 0, 

fij = — £ijk —3 0 {—xa) . 

y x 

The general properties of the steepest descent configuration that minimises the action 
in eq. ( p.7[ ) can be understood using the following argument. The insertion of the operator 
M into the path integral annihilates a magnetic monopole at time t = 0. Thus one 
expects that the leading configuration has a smooth magnetic charge density which at 
negative times integrates to a different value than at positive times, the difference being 
the magnetic charge of the monopole annihilated by M. One might expect that the action 
is minimised by a spherically symmetric configuration. The natural candidate for such a 
magnetic charge/current density would be 

3n( x ) = - 7 % > ( 3 - 9 ) 

g |xp 

since it has a divergence at the point in space-time where the monopole is created, and thus 
is naturally associated with a space time event of annihilation of a monopole. The overall 
coefficient is determined by the requirement that the 3D flux emanating from this point is 
the same as that created by M. Such a current cannot however be constructed from the 
field strength at all, since by definition any magnetic current that can be written in terms 
of the fundamental fields of the theory is divergeneseless. Thus one needs to balance the 
total flux that emanates from the point x = 0. The simplest and the most natural way to 
do that is to have the balancing flux going in an infinitely thin flux tube: 


= 
J fl 


1 


Jl - 47T <5^4 5 3 (x)9(-x 4 ) 


(3.10) 


The thin flux tube term in eq. (3.1C) has a natural interpretation of the world line of the 
monopole annihilated by the operator M. This structure is very similar to the one arising 
in the calculation of the VEV of the magnetic vortex operator in 2+1 dimensions [7] and is 
generic for a broken symmetry phase. Pictorially it represents the flux of the spontaneously 
broken charge emanating from a single space-time point in a spherically symmetric way, 
like the head of a dandelion, while the necessary “nutrients” (incoming flux) are supplied 
by a thin stem extending infinitely into the past. We will thus refer to a configuration of 
this generic type as a “dandelion”. 

We will find however that a spherically symmetric dandelion does not minimise the 
classical action. Instead the action is lower for an axially symmetric dandelion, the spherical 
symmetry being broken due to the sensitivity to the direction of the Dirac string. In 
the rest of this section we will construct the vector potential for an axially symmetric 
dandelion and consider the minimisation of the action. Although we cannot prove that 
the configurations of the type we consider lead to an absolute minimum of the action, we 
believe that parametrically (and physically) our results are correct. 

We start with the following expression for an axially symmetric vector potential and 
the Higgs field in four Euclidean dimensions: 


H a = 


4 a_/(^+>„ X b 

w lib o 


(3.11) 
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where 


2 2 , 2,2 

r — x 1 + x 2 + x 3 , 


x 4 Q 

z = — , X = cos d . 


Calculating the components of the field strength we find 

njk 


Xi df x 3 Xi df\ 


where 


/7“ — /?.. if a _l I 
r i3 ~ r% 'J n ^ gr 3 

T?a a A a 1 xb a ft ^ 1 Xk d f 

F 4i = d 4 Ai = -e aib ^d 4 f(z,x) = - e aik jr-, 

g r z g r 6 oz 




g r° 


The covariant derivative of the Higgs field is given by 

( Di H) a = -{S ai -^)(l-f). 

ryt /y* Zj 

The magnetic current for this configuration is 

1 


•M = 
J (1 


g 


TlC 1 _ _ 4?n W S 3 (x)6(-x 4 ) 


(3.12) 

(3.13) 

(3.14) 

(3.15) 

(3.16) 

(3.17) 


This looks like a radial dandelion but with the density of the current not spherically 
symmetric - a dandelion in the wind. 

With these expressions we find that the classical action in eq. (3.7) is given by the 
following integral 

r-oo dr 


with 


Si = 


Sdandelion — 


dz 


IT 

g 2 Jo 


[Si + S 2 + S 3 + S 4 + S 5 ] 


f 1 dx [ 2 / - f 2 - 0(—z) 

. 7—00 J — 1 

S 2 = Il dz L dx {^) • 


(3.18) 


(3.19) 


S3 = 


dx 2 z 2 


d£ 

dz 


/ OO f‘ 1 

dz / ( 

-00 J -1 

Si= H dz L dx(i - x2) {^) 


s 5 = 


dz f dx 0(—z)(l — x)(l — /) 

. J -1 


—-—(1 — /) — 2^- 
1 + x v cte 


The term Si comes from the square of the abelian component of chromomagnetic field, 
the term S 2 from the chromoelectric field, while S 3 + S 4 + S 5 arise due to the square of 
the nonabelian part of the chromomagnetic field, the square of the term in eq. (377) which 
involves the covariant derivative of the Higgs field, and the interference term between the 
two. Defining 

K = l-f, (3.20) 
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with A an ultraviolet cutoff, and L the linear size of the system. 

The ultraviolet divergent has to do with the fact that our monopole has a zero core 
size, and is therefore not a problem. The infrared divergence appears since our ansatz for 
the dandelion vector potential is dilatationally invariant. As we will argue shortly, this 
is unavoidable and indicates that in perturbation theory the expectation value vanishes 
if the volume of space time is infinite. We believe this is an artifact of the perturbative 
calculation, and nonperturbative confinement effects must provide an effective infrared 
cutoff on this calculation, regulating this divergence. 

The most important question is of course whether the remaining x and z integrals are 
finite. It is easy to see that this is the case for any smooth function which satisfies the 
following boundary conditions 

K(z, — 1 ) = 0 , z < 0 , 

I K(z —■> oo, x) — 1| < — , (3.22) 

z 

| K(z —► —oo, x) | < — . 

We take as an ansatz for K functions of the form 

K(z) = p(z)( 1 + x) 1 "^ , (3.23) 


where p(z) is a smooth approximation to the step function. As an example consider the 
one parameter class of functions 

p(z, 7 ) := —b — arctan (— - = ) . (3.24) 

K 2 vr Vv/1 + 7 W ^ ' 

From these we can construct a two-parameter class of functions K(z,x,a, (3) where 

K(z, x, a, (3) = p(z, a)(l + x) 1_p ^’^ . (3.25) 

7f 

Writing dandelion (ot, (3) = —tt ln(A L)c(a,(3) and calculating the integrals numerically, we 

9 

find as shown in Figure 1, a minimum value for the action at a ~ 1.35 and (3 ~ 2.65 
corresponding to c(a, (3) ~ 2.17. Several other smooth approximations to the step function 
have been investigated and all lead to a higher action. 

A striking property of these boundary conditions is that a spherically symmetric ansatz, 
corresponding to the function K independent on the angle x has an action which diverges 
worse than logarithmically. In particular if we take 
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a 



Figure 1: Plot of the two parameter variation c(a,(3) for the dandelion action. 


which gives the spherically symmetric magnetic current eq. ( |3.9| ), the z,x integral diverges 
as the second power of the infrared logarithm. The spherically symmetric configuration 
thus gives a negligible contribution to the path integral for the expectation value of M. 


Physically the meaning of the dandelion configuration is quite clear. As X 4 —> — 00 , 
the magnetic field eq. ( |3.13| ) is the field of the point-like ’tHooft-Polyakov monopole. At 
negative but non-infinite times an extra magnetic charge density appears at distance scales 
of order X 4 . This extra charge density is not spherically symmetric but rather has only 
axial symmetry and is pushed away from the Dirac string. Since our ansatz eq. ( 3.1 1| ) 
is dilatationally invariant, we expect that for any field of this type the extra magnetic 
charged density is small within the cone of the opening angle X4 Jr around the string. As 
X 4 —* 00 , the field vanishes. Thus from this perspective it looks like a configuration that 
interpolates between a monopole in distant past and a vacuum in distant future. Note that 
even so, the conservation of the magnetic current cannot be violated by the dandelion. To 
see this explicitly we have to analyse the field profile at fixed X 4 . As long as r < \x±\, 
the field is what we have just described - a monopole at negative time, and vacuum at 
positive time. However at r > \x^\, the field is given simply by the value of the function 
K at z = 0. As long as K is a continuous function of z, this means that the magnetic flux 
at spatial infinity is time independent. The balancing flux is precisely due to the extra 
magnetic charged density which lives at distances of order r oc x 4. Thus the dandelion 
configuration describes at negative times a point-like ’tHooft-Polyakov monopole, which 
at radius r ~ \x±\ is surrounded by an axially symmetric “shell” of negative magnetic 
charge density. As \x±\ —> 0 the radius of the shell shrinks to zero. At positive times the 
configuration has no core but instead all the magnetic charge density is concentrated in a 
shell whose size grows as X 4 


00 . 


Note that the logarithmic infrared divergence is not an artifact of our dilatationally 
invariant ansatz. It is easy to see that one cannot do better than that. Let us relax the 
assumption that K depends only on the ration X 4 /r. A little algebra then leads to the 
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following expression for the action of the dandelion configuration 


'S'dandelion — ^ J j 4 j dx 

- x dr dx 4 / dx 

9 J J - 00 J -1 


1 . r , 2 , 2 fdK\ 2 (dK 

3(1 -Vi 2 +( —) + 2 (— 


9 x 4 y 


1 






OK \ 2 / 0 KV 
0 X 4 y \ 0r y 


(3.27) 


11-x/ 0iY 

+ ^TTx ( A + (1 + x) 0x 


Let us follow the change of K with time X 4 < 0 at fixed value of r. The qualitative 
behaviour is clear from our discussion above. At large negative times K = 0, while at 
times close to zero, K is a fixed function k(x) which does not vanish at large r. Let’s call 
the time interval during which this change in behaviour occurs A T(r). We can estimate 


the magnitude of various terms in eq. (3.27). In particular 


1 4 A T(r) 4 

ax^—^K oc-7T —k , 


(3.28) 


’ —00 

f 0 


dK 


1 




V 


It is obvious from these two expressions that T(r ) —^r^oo Ar, or else the subsequent integral 
over r diverges stronger than logarithmically. 

Thus we conclude that in the leading order in perturbation theory the expectation 
value of the monopole operator is 


7TC . 


(M) = exp (- - ln(AL) 


(3.29) 


with constant c ~ 2.17. 


4. Discussion 

In this paper we have shown how to set up the perturbative calculation of the expectation 
value of the monopole creation operator in the confining phase. Perturbatively we find that 
the VEV vanishes in the infinite volume. However our result eq. ( |3.29| ) is an intrinsically 
perturbative one, and as such we certainly do not expect it to stand beyond perturbation 
theory. Nonperturbative contributions will provide an infrared cutoff. It is thus very likely 
that the (M) 7 ^ 0 in the confining phase. 

A perturbative calculation of VEV is of course not reliable since it is dominated by 
the infrared physics. A quantity which is more amenable to such a calculation is short 
distance behaviour of the monopole-antimonopole correlation function [12]. Although we 
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have not calculated this quantity, the perturbative cutoff dependence of (M) suggests that 
the correlator behaves as 

(^M{x)M^ (y)^ oc |x — y\ s* . (4.1) 

This indicates that even in the pure Yang-Mills limit the monopole operator does not freeze, 
even though its phase must decouple from the finite mass spectrum. Instead M couples to 
finite mass excitations via its modulus. 

It is worth noting that the present approach can be easily extended to finite temper¬ 
ature. It is interesting to see whether at high temperature the VEV of M vanishes in 
perturbation theory. 
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